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In this document, an attempt to derive a cardinalized probability hypothesis density filter 
(CPHD) (corrector formula only) for extended targets is given. 

1 Derivation 

We first give some quantities related to the derivation. 

• The extended target measurements are distributed according to an iid cluster process. The 
corresponding likelihood is given as 

f(Z\x) = n\P z (n\x)Y[ Pz (z\x) (1) 

zez 

• The false alarms are distributed according to an iid cluster process also. 

f(Z FA ) = n\P FA {n) J] Pfa(z) (2) 
zez FA 

• The multitarget prior f(X k \ZQ :k -i) is assumed to be an iid cluster process. 

f(X k \Z 0:k -i) = n\P k+l \ k (n) Yl Pk+i\k(xk) (3) 

x k <=X k 

where 

Pk+i\k(xk) = N^l^Dk^ixk) (4) 

with TVfcifc.! = /I>fc|fc_i(xfc) dx k . 
The p.g.fl corresponding to the updated multitarget density f(X k \Z k ) is then given as 

G Mk [h] = 5Z x k (5) 

where 

F[g, h] 4 J h x G\g\X\f{X\Z k _{)8X (6) 

G[g\X] 4 f g z f(Z\X)5Z (7) 

• Calculation of G[g\X\: Suppose given the target states X the measurement sets corre- 
sponding to different targets are independent. Assume that targets are detected with 
probabilities pd(x). Then the p.g.fl for a single targets measurements become 

G[g\x] = l-p D (x)+ PD (x)Gz[g\x] (8) 

Then p.g.fl for the measurements belonging to all targets become 

G[g\X] = (l-p D (x)+p D (x)Gz[g\x]) x (9) 

With the addition of false alarms we have 

G[g\X] = G FA [g](l-p D (x)+p D (x)G z [g\x]) x (10) 
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• Calculation of F[g,h]: substituting G[<?|X] above into the definition of F[g, h] we get 

F[g,h] = J h x G FA [g](l-p D (x)+p D (x)G z [g\x]) x f(X\Z k _ 1 )SX (11) 

=G FA [g] J {h(l - p D (x) + p D (x)G z [g\x])) x f{X\Z k ^)5X (12) 

=G FA [g]G klk _ 1 [h(l-p D (x)+p D (x)Gz[g\x]] (13) 
=G FA [g]G klk _ 1 [h(l-p D +p D G z [g}] (14) 

where we omitted the arguments x of the functions. 
For the iid cluster processes, we know the following identities 

G FA [g] =G FA {p FA [g\) (15) 
Gz[g\x]=G z (p z [g}\x) (16) 

where the functions G FA (-) and Gz{-\x) are the probability generating functions for the car- 
dinality distributions P FA {-) and P z (-\x) respectively. We also know the following derivative 
expressions. 

±G[g] = GW(p[g])l[p(z) (17) 

zez 

Now taking the derivatives of F[g, h] with respect to Z. 

—F[g,h] =G { F A (p FA [g})p FA (zi)G k \ k _ 1 (p k \ k _ 1 [h(l - p D + p D Gz(p z [g]))] ) 

+ G , F J 4(PFA[5])^J._ 1 ^fe|fc_i[/i(l -pd +PDG z {Pz[g]))]^Pk\k-i[hpDG { z\pz[g])Pz{zi)] 



=F[9MpFA{zi \ G F A{p FA [g]) 



G k \ k -i{Pk\k-i[K l - pd +PDG Z (Pz[g}))] 



+ 7 : -\Pk\k-i hp D Gy(p z [g\) 

G^-^p^iKl-pD+pDGzipzig}))]) 1 P^( z i) 

=F[g,h] PFA (z 1 )^ A (p FA [g]) 
+ ^k\k-i{Pk\k-i[K l ~ pd + PDG z {Pz[g]))])pk\k-i hpDG{ z\pM) ^^^ 

where we use the functions Cfa(') and Cfc]jfc_i(') defined as 



(18) 



^ A (§M)" 1)=,OS " ,aFAx) <I9) 

(1) , . \ 
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The second order derivative is given as 



^-F[g,h]=F[g,h]( ]J p FA {z')\ ( ]J U^faM) 

Z'=Z\,Z2 \z'=Zl,Z2 

+ Cfc|fc_ X (pfc|fc-l [M 1 ~PD + PDG Z (Pz[g]))]^jPk\k-l 

+ Cf A {PFA[g]) + C { k 2 l_ 1 {p k \ k -i[h{i -pd +pdG z {pM))]) 

x II Pk\k-i[hp D Gfij)M)^X 

z'=z 1 ,z 2 

+cl}L(^| fc -iKi-^+p D G Z (p z [ 5 ]))]) Pfc | S; _ 1 [^4 2) feb]) n 

yl — •r-i Vr> 



(21) 



Then we can write the general formula as 

±F[g,h]=F[g,h}(Hp FA (z'))Y^ U feWbD 

+ E ^(Pfclfc-i^l -PD +PDGz{Pz[g]))]) II Wl^l) (22) 
QZW VeQ ' 

where 

=Pk\k-x[hpDGf U(p g \g]) ]J ^^y] (23) 

A proof of ()22p can be found in Section [6l Setting 5 = and taking derivative with respect to 
x. 

Jx~Jz m k] = Jz m h]C ^-l i^- 1 ^ (1 ~ PD+ P dGz ^ ) 
x (1 -p D (x) +PD(x)G z (Q))p k \ h _ 1 (x) 

+F[o,h](i[p FA v))Y,( II (cS^co) 

^z>ez ' vzz \ w&v ^ 

+ e ^(pkifc-itMi-PB+p^zco))]) n w[o,^) ) 

I Eqzw Cm-i (Pk\k-i [Kl -PD+ PdG z (0))] ) UveQ W [0, h] 



E 

WdV 



Cff 3 (o) 

+ Egzw Cl'i^ (pfc|fc-i [*(1 -PD+ pdG z (P))] ) IlyeQ W [0, h 



(24) 



: — F[0,fc]cgL(p*l*-i[^ -pd +p D G z (0)) 
x (l -p D (x) +PD{x)G z (0))pk\k-i( x ) 

+F[o, h ](i[p FA ^))j:( n (^a d (o) 

+ E cf|S(p fc | fe -i[Mi-PD + PDGz(o))]) n V v[0,h]\ 

QZW VeQ ' . 



/ E Q z^(ny eQ w[0,/ i ])(4jS 1) (p fe |^i[Ml-PD+PDG z (0))^ N 
(1 -p D (x) + PD{x)G z (0))p k \k-i(x 



+df2APk\k-i[h(l-PD+PDG z (0)) 



\ 



x l^VeQ vv[0,h] 



J 



+ Eqzw CkS-i (Pfclfc-i [^(1 " Pd + PdGz(O))] ) ny eS W [0, /»] 



(25) 



Evaluating at h = 1, we get 
5 5 (5 



— — F[0, 1] =^[0, l]C$_i (pfc|fc-i [1 - PD + PdG z (0)} ) (1 - Pd{x) + PD{x)Gz{0))p k \ k - 1 {. 

^'e2 ' -pzz \wev ^ 

+ E C^^-ifl-PD+PDGzW]) [J W[0,1]) ) 
qzw vgS ' / 



x (1 - p D (x) + pu(a;)Gz(0))p fc | fc _i(a;) 



xE- 

WGP 



X 2-VgQ »jv[0,1] 



C^ l} (0) 

+ Eqzw Ck\?-i {Pk\k-i [1-PD+ PdG z (0)] ) UveQ W [0, 1] 



Defining additional quantities 



(26) 

(27) 
(28) 



QZW 



we get 

^^F[0, 1] = ^[0, lRifL (pfc|fc-i [1-PD+ PdG z (0)] ) (1 - p D (x) + p D (x)G z (0))p klk _ 1 (x) 

+ Fp, i] ( n pfaw) e ( n M 

( T,Qzw a Q( y (k\k l -i ) (Pk\k-i[^-PD +PdG z {0)] 

x(l-p D (x) + PD(x)Gz(0))p k \k-i{x) 
(3 W - PD +PdG z (0)]^ ^ 29 ) 



E 



\ 



V 



X ^VeQ ri v [0,l] 



J 
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Now dividing by jgF[0, 1] to obtain D k i k (x), we get 

D k \k(x) =(1 ~Pd(x) +PD(x)G z (0))(^l_ 1 (p k \ k _ 1 [l -p D +PDGz{$)])p k \k-i{x 

( Eqzw a Q (CkS-^ (Pk\k-i [l ~~ Pd + PD 

Y.vzz [ Tlwev Pw ) J2wev ^ 



x (l - p D (x) + PD{x)G z {tij)Pk\k~i{x) 

k\k- 



H ( k \flAPk\k-l[^-PD+PD 



+ 



Defining coefficients oj-p as 



x 2^veQ v v [o,i] 



Y^vzz n 



Pw 



UJ-p 



n 



Wev 



Pw 



J2vzz Tlwev 



Pw 



(30) 



(31) 



we get 



D k \ k (x) ={l- PD {x) +p D {x)G z {0))C { k 1 l _i (pk\k-i [l ~~ Pd + PdG^(O)] ^)p k \ k -i(x) 

' T. Q zw^Q{c^ l \pk\k-A l -PD+PDG z ^)]) N 



■pzz VKeP 



/3w 



x(l -pd(x) + P.o(x)Gz(0))p jfc | fe _ 1 (x 
>fc|fc 



+C { k \?l 1 (Pk\k-l^-PD+PDGz(0) 



X 2-VGQ ^Joo] 



(32) 



= (1 -p D {x) +PD(x)G z (0))(V_ 1 (p klk 

' Y.QZwUQ(c^X l \PHk-l[l-PD+PDG Z (V)]) ^ 



■pzz vkgp 



x(l- PD (x)+p D (x)G z (0)) 
+CfciS- ) i (Pfc|fc-i [l — Pd + PdG z (0) 



»7v[0,l] 



If now we define a constant 



K -E^E^E aeCfcjSJ 15 (Pfc|fc-i [1 +Pr»Gz(0) 



n/z ireP ^ ©zn 



Pk|fc-i(af) 



(33) 



(34) 



we get 



Afc|fcO»0 =(Cfcffc_i + «)(! -Pafa) +PD(a;)G z (0))p fc | fe _i(x) 



»fc|fc- 

-pzz wep ^ W ©zw ve© 



p D (x)G { P(0) p 2 (z'|x) 



n 



^[0,1] j^v PFA ^ z 



jTPk\k-i(x) (35) 



where we have dropped the argument of the terms Cjt|fe-i {lPk\k-i [l ~~ Pr> +PdG^(0)]^ for sim- 
plicity. 
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We have G k \ k [h] as 



Gk\k[h] 



\ 



+ Eqzw ( { k \k-i (Pk\k-i [Kl -PD+ PdG z (0))] ) UveQ Vv [0, h] 



\ 
J 



F[0,l]EvzzUweV (C?T D (0) 
T,qzw Ck\?-i {Pk\k-i [1-PD+ PdG z (0)] ) UveQ Vv[0, 1] 
G k \k-i (l>k\k-i [K l -PD+ PdG z (0))] ) 



(36) 



))]) 

G k \k-i [Pk\k-i [1-PD+ PdG z (0)] ) 
ZvzzUwev Ufa 1 } (0) + E Q zwCk\k-i{Pk\k-i[h(l - Pd +PdG z (0))} ) UveQ w[0, h] 



J2rzz Uwev 



(37) 



Then, 



Gfc|fc(a 



Cfcifc-i (xp k \k-i [l-p D + PdG z (0)] 
G k \k~i (Pk\k-i [1-PD+ PdG z (0)] ) 

YZvzzUweV UfT\°) + Y^qzw a Q xlQl C ( k\?\( x Pk\k~i[ 1 ~ Pd + PdG z (0) 



Y2vzz U 



vzz iiwer 



(38) 



G fc | fc _i (xp^fc.! [ • ] ) Epzz IW (47(0) + E S z^ ^'^CJiS?! (*P*|*-i [ 



G 



k\k-l 



Pk\k-l[-]) 



(39) 



where we skipped the arguments of Pk\k-i — Pd + PdGz(0)] for the sake of clarity. 
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2 Summary 

Calculate the quantities D k \ k (x) and G k \ k (x) from the given quantities D k \ k _i(x) and G k \ k ^\(x) 
as follows. 

D k \ k (x) =(4|fe_i + K)(l -p D (x) + PD{x)G z (0))p k \ k _ 1 {x) 



G k \ k (i 



+ e «p e £ e ^ e yD( ;!'y 11 (0) n < 4 »> 

G k \k-i(xPk\k-i[-]) ^vzzWw^v U { FT\°) + T, Q zw a Q xlQl C { k\?\( x Pk\k-i[ 



Gk\k-i{pk\k-i[ • ]) 



Z~2vZZ Y[\V£V Pw 



(41) 

where we skipped the arguments of p k \ k _i [l — Pd + PdGz{0)] for the sake of clarity and 

Pk\k-i(x) = n fc|fc ~ 1 ^ (42) 



« - E E J" E "eCi'iS" (»!*-' I 1 - m + pdG z (0)] ) (45) 

VZZ WeV PW QZW 

UV ff R~ ( Ab > 

l^vzz l lw<=p 

Pw =4V ] (V) + £ CkS-i{pk\k-i[l-PD +PDG z (0)])a Q (47) 

QZW 

veQ 

Wl9, h] % fc | fc _! [/ipsGf'W]) n ^7j] ( 49 ) 
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3 Calculating Posterior Cardinality Distribution Pk\k{n) 

In this section, we are going to calculate the posterior cardinality distribution Pk\k{ n ) based on 
the posterior cardinality probability generating function (p.g.f.) G^[x). For any cardinality 
distribution, p.g.f. pair given as P(n)-G(x) we know the relationship 

P(n) = i-G^(O) (50) 

TV. 

which we are going to use in this calculation. In the previous sections, we calculated the posterior 
cardinality p.g.f. G^ix) as 

G*| fc -i (xp*| fc -i [ • ] ) Zr*z Kwev (cff ) (0) + J2 Q zw «s^ |S| C?|£ [ ■ ] ) 

^ ) G fc | fc _ 1 (p fc | fc _ 1 [-]) EvzzUwevPw 

(51) 

Here, the denominators on the right hand side are constant with respect to x and for this reason, 
we are going to take the derivative of only the numerator terms. 
For the derivative of the multiplication of two functions, we have 

n 

{f(x)g(x)) in) = Y J C n ,f {n ~ l \x)g {%) {x) (52) 

i=0 

where C n ^ is the binomial coefficient given as 

c n\i = i]{n i l)v ( 53 ) 

We have the generalization of this derivative rule as 

M \ (n) M 



UfiW) = E C nlil _ iM Hf^\x) (54) 

K j=l J 0<ii,i 2 ,...,i M <M j=l 
il+i2-\ ViM=n 



where C n ^ 1>mmm> i M is the multinomial coefficient given as 



C n \h,...,i M — ■ ■ • , . . ■ (55) 
We are going to use the following lemmas in the calculation. 
Lemma 1 Derivatives of the term {xpk\k-i [ ' are given as 

G%-± (*P*|*-i [ ■ ] ) L =0 = [ • ] ) "^-i(0) (56) 

=n\(p k \k-i[-]) Pk\k-i(n) (57) 

where Pk\k-i(') i s the predicted cardinality distribution. □ 



Lemma 2 



Derivatives of the term ((p^\o) + YIqzw a Q x ^(j!\k-i ( x Pk\k-i [ ■ ]) J are given 



as 



4V\V+ E ^ ISI CSa(^| fe -i[-]) =C?T l) (0)<o + n! £ a c <g2> l( 0) 

QZW J x=0 QZW 

\Q\=n 

(58) 



Notice that the derivative is equal to if n > \W\. □ 
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Lemma 3 The nth derivative of the summation 



e n (c^ i} (o)+ e ^^'ciia^-ic-])) (59) 



is given as 

\V 



E E Cnln,...,^, II (0)^,0 + E «flCi|ffi(0)) (60) 



VHWj} 1 ^ i 1+ i 2 +...+i ]Pl = n 



VZZ ^ 0<ii,i 2 ,...,i|p|<|P| J=l 
n! 

•PZZ 0<i 1) i 2 ,... ) i|- P |<|7'| 

^TX i-, J-i»J 1 



E E jn[(cS? l) (o)^,o+ E «scil?i(o)) (6i) 

PZZ 0<ii,i 2 ,...,i|p|<|P| i=i V S z ^j ; 



VHWj}fX i 1+ i 2+ ... + i m = n \3\=ij 



i,<\Wj\ Vj 

where we used and the condition i,- < \Wj\ Mj of the second summation is added just to 
reduce the number of possibilities. □ 
Now we are ready for the final result using the derivative of the multiplication rule ()52[) and 
the three lemmas above as follows. 



\ n n—i 

Pk\k(n)= j — r ^2(Pk\k-i[-]) Pk\k-l( 

Gk\k-1 \Pk\k-l [ ' J ) YjVAZ Ww&V @ w i= 



n — i) 



E E nfe'"*..+ E (^) 



ij<\Wj\ Vj 

4 Reducing CPHD for extended targets to PHD for extended 
targets 

Now we check whether the CPHD equations given above reduce to PHD when the cluster 
processes are replaced by Poisson processes. This is done also for checking the formulas given 
above. First we give the related quantities. 

• The extended target measurements are distributed according to a Poisson process. The 
corresponding likelihood is given as 

f(Z\x) = n\P z {n\x) \\ p z (z\x) (63) 



where 



P g (n\x) = ~ (64) 



The false alarms are distributed according to a Poisson process also. 

f(Z FA ) = n\P FA {n) J] p FA {z) (65) 
zez FA 



where 



Pfa(u) = P- (66) 

n! 
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• The multitarget prior f(X k \Z . k _i) is assumed to be a Poisson process. 

f(Xk\Zo:k-i) = n\P k+1 \ k (n) Y[ Pk+i\k(xk) (67) 

where 

Pk+i\k{xk) = Nj^l^D^^Xk) (68) 
with N k \ k _i = f D k \ k _i(x k ) dx k . We have the cardinality distribution P k +i\ k (n) given as 

^fc+i|fc(n) = ij^ 1 (69) 

We here must note that the probability generating function and functionals corresponding to a 
Poisson process with parameter A and density p{) corresponds to 

G(x) = exp(Ax - A) and G[h] = exp(\p[h] - A) (70) 

We have 

(x) = X n G(x) = A n exp(Ax - A) (71) 
Also the following expression holds. 

G (n) (0) = X n e~ x (72) 
The functions C,fa and Cfelfc— l are then given as 



> = (a)(*- i )= r ' ; A^., (73) 



i > 1 



\^Gfa(^) / 1 A i = 1 

Cff*_iN = ( F^TTrl = = Jv ' > !=^-Ai (74) 

■"' 1 Gk\k-i{x)) / I JVfclfc-i * = 1 



Then in the formula 
_6_ 

Jz 



F[g,h]=F\g,h](l[p FA (z'))Yl U (cST'WbD 

+ e cfiS^-itMi-PD+PDGz^b]))]) n^^)' 

qzw veQ ' 



(75) 



all terms in the summation X]qzw becomes zero except Q = {W} which means that there is 

only one set in the partition Q and it is V = W . Also from above, we have CjJa' ( x ) = ^|WI,i- 
Substituting these into the above equation 

F[g, h] =F[g, h] ( J] Pfa(z')) £ II ( A Vl,i + N k\k-iVw[g, h]\ (76) 

This equation is the same as the equation Mahler derived for 4^F\g, h] [1, Eq. (28)]. In the 
CPHD equation 

D k \ k {x) =(Cfc}fe_i + ^H 1 ~Pd(x) + PD(x)Gz{0))p k \ k -i(x) 

+ \ <*> L fa L L w[0 ,i] II z^y^-i(-) (77) 
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By above facts, it is easy to see that n = 0. Substituting above facts into this equation we have 



D k \k(x) =(l ~Pd{x) +PD{x)G z {ti))D k \ k _ l (x) 

^ \- U{W} pp(x)G^ Wl) (0) 1 T p z (z'\x) 



■pzz ive? 
Seeing that ctryy\ = rjw[0, 1] we get 



(78) 



z'€W 



D k \k(x) = (l ~Pd{x) + PD(x)G z (0))D klk _ 1 (x) 

+ 2^ ^ 2^ 11 zzj^ d ^-^ 



VIZ WdV 



z'&W 



Pfa(z' 



We also see that 



Pw =C^ I) (0)+ $*-i{Pk\k-i[l-PD+PDGz{0)])a Q 

QZ.W 

=^\w\,i + Nk\k-mw[0, 1] 



=M\w\,i + N k\k-iPk\k-i pdG^z \°) 1 [ 



Pfa(z') 



(79) 



(80) 

(81) 
(82) 

(83) 

(84) 
(85) 



Dividing and multiplying all quantities in the last summation by \\ w \ and distributing, we 
obtain 



D k \ k {x) =(l - p D (x) + p D (x)G z (0))D klk _ 1 (x) 

pp(x)G^ Wl) (0) -r-r p z {z'\x) 



•pzz w&v 

Defining the new coefficient 

A „ 



(86) 



(87) 



^ w+ %- I K l \o)n^| ) 
=vu+^ifc-i^4 wi) (o) n 



Apfa(^) 



(89) 



which are the same coefficients defined in pQ. Also it is easy to see that 



UJ-p 



J2vzz Ylw&v Pw 



J2rzz Tlwcv 



dw 



(90) 
(91) 
(92) 



12 



which are the same u-p coefficients in pQ. Knowing that G^'(O) = j^ w \(x)e we get 



D k \ k (x) =(l- PD (x) + PD (x)e 7(x) )-D fc | fc _i(x) 

+ 2^ 2^ ^ 11 



which is the same formula in [1, Eq. (5)]. 



z'ew 



Xp FA {z 



(93) 



5 Reducing CPHD for extended targets to CPHD for standard 
targets 

For standard targets, P z {n\x) = except for n = 1 where it is unity i.e., P 2 (l|x) = 1 which 
makes Gz{x) = x. Then 



W[0,h] ± Pk \ k -i\hp D G(P(0) [J -^p- 
L Xv Pfa{z) 



'Pk\k-i[hPD^j\ \V\ = l,V = {z'} 
otherwise 

Pz(z') 



-Pk\k-1 



hpD 



Pfa(z') 



>W\,1 



This gives 



«s - n w[o, i] = jj Pk\k~ 

V&Q z&W 



PD 



Pz(Z ) 

Pfa(z') 



\Q\,\W\ 



where Q is a partition of the set W. 



(3w =C?T l) (0)+ Yl tkS l \(pk\k-i[l-PD+PDG z (0)})a Q 



QZW 

<;?(»)+C(ftM[i-PD]) n^-i 

z'ew 



PD 



Pz(z') 



Pfa{z') 



(94) 

(95) 
(96) 

(97) 

(98) 
(99) 



Substituting these into GW[/t] we get 
A Gk\k-l (Pk\k-i [h(l ~PD + 

Gk\k[h] 



G k \k-\ (Pk\k-i [l ~~ Pd + PdG z {0)} ) 

C£|S (Pk\k-i [h(l -PD+ PdG z (0))] ) IlygQ W [0, , 



Gk\k-i\Pk\k-i[K l - Pd) 

Gtyk-l \Pk\k-l [l ~PD 

Evzz Uwev ( (fa + Cfc'K (Pk\k-i [h(l - Pd)] ) U z 'ew Pk\k- 



(100) 



hp D ^kX 

^ u PfaKz') 



Evzz Ylwev Uf7\°) +Cfcpi(p*|ft-i[l-P£>]) II*'ewP*|fc- 



(101) 
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G k \k-i\Pk\k-i[K l -Pd)] 
Gk\k-l[Pk\k-l[l - Pd 





logG FA (p FA [g}) + log G fe |fc-i (pk\k-i [h(l -PD+PDPz[g})]^ 


) 


9=0 




l°g GfA {PFA [9] ) +logG/ c | fe _i(p A; | fc _i[l -PD +PDPz[g]]^ j 




5=0 



(102) 

G k \k-\[Pk\k-i[K l ~Pd) 



G k \k-i(pk\k-i[l -P£>]) 

EpZZ UweV 5W lo § ( GfA (j?FA [d] ) 1 (Pfc|fc-1 [M 1 - + PDPzbD] 



3=0 



T,VZzYlw€VSW lo & [ GFA(pFA\g])G k \ k _ 1 (p h \ k _ 1 [l - p D + PDPz[g\] 



9=0 

(103) 

where we used the facts that J (p A = logGFA and J 1 = l°g 1 to obtain (|102p . We 

can see that for any p.g.fl., G[g] 



e n (^M+^irf) 



from the product rule for functional derivatives (See [2 8th row of Table 11.2].). Since G[g] is 
arbitrary, we can replace it with G[g] — G[0] which gives 

£ n ((obi - ^r 1 + §a) - Sim -• S[01)IZ ' 005) 

Now, evaluating both sides at <? = 0, we get 



(106) 

9=0 



Now suppose we have two functionals and we need to evaluate the following ration. 



8G1 r«n g(Gi[g]-Gi[0])' z ' 



^ (107) 



This division gives jj indeterminate form. One must use l'Hopital's rule \Z\ — 1 times (take 
derivatives of the numerator and denominator \Z\ — 1 times with respect to g), to get 



J2vzz Tiwev & sw [2] 
S-pzz Eli^eP §jy[0] 





9=0 


<5Z 


<5G 2 [<?] 


9=0 


<5Z 



(108) 
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Using this formula for evaluating (|ll)3p . we obtain 
G k \k-i{vk\k-i[K l ~Pd)] 



Gk\k[h] 



Gk\k-l[Pk\k-l[l ~PD 



f^log ( G F A(PFA[g])G k \ k _i(p k \ k _i[h(l - pd + PDPz[g])] 

x ± -^-!=* ,109) 

JZ {GFA\g]G k \ k -i[h(l-p D +PDPz\g\)]) 



9=0 



9=0 



jz [GFA[g}G k \ k -x[i-PD+PDPz\g])) 



(110) 



9=0 



which is the formula for CPHD for standard targets (See [5J Equations (113) and (114)].). 



6 Proof of the Main Equation (E2D 



The proof is by induction as in pQ. The formula can be seen to be satisfied for first and second 
order derivatives in (fT8j) and (f2T|) respectively. We assume that for \Z\ = m the formula is 
satisfied as below. 

^F[g,h] =F[g,h]( ]J p FA (z')\ £ J] U F ^\p FA [g]) + £ Cj^i II ^I^l) ( U1 ) 
Now we need to show that the formula is satisfied for Z U z+. 



SzZSZ F[9 ' h] =1L 



F[ g ,h]( n ^a(.o) e n tow^D + e cgia n wi^i) 

(112) 

*m( II p^(«0) (tikpFAlg]) 



+ Ci[jfc_i»/K}bi'»] 



*e n te'^bD+ECn^^i) 

p/zwe? v qzw veQ ' 

+F[ 9 ,h]( n PFA (/)) e f n fe i WbD+ e ciis?! n 

^ \ +£ Q z W (nv eQ wM)^ 
iy G p c?T l} (PFAb]) + Eqzw ci'S-i riv-Gs wb, fc] 

(113) 

-F[ 9 ,h]( n p ^(/))e n (csr^bD+E^iiiwb^] 

^-z'ezuz* ' v zz w &v\j{zS\ ^ szw veQ 

+i^]( n ^^)) e ( n U F A i \pFA[g])+ e ci',s?i n 

^z'ezuz* ^ pzz VweP ^ qzw VeQ 
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+ E e z l ,(n, eQ w[^])E, eg ! ^f 1 ) / 

^ A\w\), r n , A\Q\) rr r '- 114 '' 



=*W n p^o)e n (cS!rv^])+Etfi?-iii^. fc ] 
+F\g,h\( n «™co) e f n (c?T i) (^b])+ e ciis?! n wi^i 

^z'ezuz* ^ pzz V we? ^ qzw veQ 



j2 ^FA l+1) (PFA[g}) + T,Qzwu{z,}Ck\?-iUveQVvi9,h} ^ 

WGP C?A D (PFAb]) + E Q ZW C { k\k-! Uv€Q WlS, /»] 



^z'ezuz* ' pzzu{z*} weP ^ szw veQ 

(116) 

which completes the proof. 
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